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We consider the vacuum energy of the electromagnetic field interacting with a spherical plasma 
shell together with a model for the classical motion of the shell. We calculate the heat kernel 
coefficients, especially that for the TM mode, and carry out the renormalization by redefining the 
parameters of the classical model. It turns out that this is possible and results in a model, which in 
the limit of the plasma shell becoming an ideal conductor reproduces the vacuum energy found by 
Boyer in 1968. 

I. INTRODUCTION 

The present paper is devoted to the discussion of the renormalization of vacuum energy in the 
presence of boundaries or singular background fields in application to the Casimir effect and it is 
aimed to partially fill the gap between the two well understood situations. These arc, on the one 
side, the Casimir force between distinct objects which is always finite and, on the other side, the 
vacuum energy in smooth background fields which can be renormalized by standard methods of 
quantum field theory. In between these two, the situation is not finally settled. Especially in [12] 
it was questioned whether boundaries can be incorporated at all into a well posed renormalization 
program. For instance, it was argued that the process of making the background field concentrated 
on a surface is not physical. 

The aim of the present paper is to discuss an example of a background field concentrated on a 
surface having both, a well posed renormalization procedure for the vacuum energy and a mean- 
ingful physical interpretation. As model we take a spherical plasma shell interacting with the 
electromagnetic field and we allow for a classical vibrational motion of the shell. The investigation 
of the plasma shell model was pioneered by Barton [3] and it is aimed to describe the 7r-electrons 
in a CgQ-niolecule. 

The heat kernel coefficients for such system are known to a large extend. Since the polarizations 
for the electromagnetic field separate into the usual TE and TM modes, one is faced with two 
scalar problems, where, however, the s-wave contribution must be dropped. For the TE modes 
it is a delta function potential on the shell. The corresponding heat kernel coefficients (including 
the s-wave) were first calculated in [7], later generalized in [10], and the finite part of the vacuum 
energy was calculated in [16]. For the TM mode the scalar problem corresponds to a ^'-potential 
on the shell and the corresponding heat kernel coefficients were calculated for a plane shell only, 
[9] . The problem with the TM mode is that the corresponding spectral problem is not elliptic and 
that the standard methods do not work. So, for example, for the plane shell even the zeta function 
cannot be defined [9]. For the spherical shell, the zeta function exists, but, as we will see below, it 
has double poles. It should be mentioned that the (^'-potential was considered in [15] (where also 
the relevant literature was collected) , however with a coupling different from that following within 
the plasma shell model (compare the Jost function in (4.41) in [15] with (18) below). 

In the present paper we calculate the heat kernel coefficients for the plasma shell model, especially 
that for the TM modes. Using these, and the simplest possible model for a classical motion of 
the shell, we construct a consistent scheme for the renormalization. Within this scheme we define 
the renormalized vacuum energy of the electromagnetic field and calculate it numerically. Also we 
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discuss several limiting cases including the limit of the plasma shell becoming an ideal conducting 
sphere. 

Throughout the paper we use units with h = c= 1. 



II. THE PLASMA SHELL MODEL AND ITS RENORMALIZATION 



We consider the plasma shell model investigated, for example, in [3] which is aimed to model 
the TT-electrons in a Ceo-molecule. These electrons are described by an electrically charged fluid 
whose motion is confined the shell. Further, the model contains an immobile, overall electrically 
neutralizing background aimed to describe the carbon atoms and the remaining electrons. The 
fluid is allowed a non-rclativistic motion. Of course, this model is a quite crude simplification, 
especially because the motion of the electrons should rather follow a relativistic dispersion relation 
[13, 14]. On the other hand side it appears to be physically meaningful and should therefore 
result in physically meaningful results for the vacuum energy. For instance, it should allow for a 
treatment of the vacuum fluctuation of the electromagnetic field coupled to the plasma shell. 

The interaction of the plasma shell with the electromagnetic field results in matching conditions 
on the electromagnetic field across the shell as shown in [3] (and earlier, for a plane sheet, in [2]). 
These conditions do not depend on the state of the excitations of the fluid. The vacuum energy 
can be calculated from the fluctuations of the electromagnetic field whereas the fluctuations of the 
fluid must not be taken into account as shown in [6] (or vice verse). In this setup, the polarizations 
of the electromagnetic field separate into TE- and TM-modes. For the electric and the magnetic 
fields the corresponding mode expansions read 

E™(t,r) = ^ 1 [e— */;,r„(fc,r)L^yj,„(i?,^)+c.c.V 

^ Jo s/2(j V VL^ J 

i™i<i 

B™(t,r) = --^ VxE™(t,r), (1) 

where L is the orbital momentum operator and uj = k follows from the wave equation . The radial 
wave function fi^m{kr) must be regular in the origin and across the shell it must fulfill the matching 
conditions 



lim (r/(,„(fc,r))'- lim (r/j,„(/c, r))' = nRfi^^{kR), (2) 

r—^R+O r— >K— 



where only the parameter 



^ 47rne^ 

^ = TTT^ (3) 



carries information on the properties of the fluid like its density n and mass m. It can be interpreted 
as a kind of plasma frequency in parallel to the plasma frequency of a dielectric. For Cgo the 
corresponding wave lengths is of the order of micrometers. The mode expansions for the TM 
polarization read by duality 

B™(i,r) = J2 r-^fe--*5i,m(fc,r-)L-i=F;,„(t?,<^) + c.c.V 
Jo 71" V2u) V VL^ J 



i>i 



E™(t,r) = VxB™(t,r). (4) 

The matching conditions are different, 

lim (rs'!,TO(fc,r))'- (rgi^m{^,r))' = 0, 

1 — >R+0 r— >K— 

lim gi^rn{kr) - lim gi,ni{kr) = ---p-(i?g;,m(A;, i?))'. (5) 

r^R+O r^R-0 k^R 
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Considered as a scalar problem, the matching conditions (2) of the TE mode are equivalent to a 
delta function potential Cl6(r — R) in the wave equation and the conditions (5) of the TM mode 
loosely speaking correspond to the derivative of a delta function. A difference is that in the scalar 
problems the zeroth orbital momentum, / = 0, or s-wave contribution is present whereas in the 
electromagnetic case it is absent, i.e., the sums over / in (1) and in (4) start from I = 1. In the limit 
fi ^ oo which is formally the ideal conductor limit the boundary conditions (2) and (5) became 
Dirichlet boundary conditions for TE polarization and Neumann for TM polarization. 

We extend this model by allowing for radial vibrations (breathing mode) of the plasma shell. In 
Cgo these arc determined by the elastic forces acting between the carbon atoms. Without going 
here in any detail we describe these vibrations phenomenologically by a Hamilton function 

2 ^ 

-f^class = 1" TT '^b (-^ ~ -^o) + Ej-est (6) 

2m 2 

with a momentum p = mR. Here m is the mass of the shell, is the frequency of the breathing 
mode, Ro is the radius at rest and £^rest is the energy which is required to bring the pieces of the 
shell apart, i.e., it is some kind of ionization energy. 

Now we consider a system consisting of the classical motion of the shell as described by -ffciass 
and the vacuum energy Evac of the electromagnetic field interacting with the shell by means of the 
matching conditions (2) and (5). We assume the classical motion adiabatically slow such that the 
vacuum energy can be taken as a function of the mountainous radius of the shell, E^^c = E^adR), 
and we neglect the backreaction of the electromagnetic field on the shell. Under these assumptions 
the energy of the classical system, £'ciass(-R) = -^fciass) and the vacuum energy add up to the total 
energy of the considered system, 

Etot = Eclass{R) + £'vac(-R)- (7) 

Next we consider the ultraviolet divergences of the vacuum energy. These are given in general 
terms by the heat kernel coefficients a„ (we use the notations of [8] and we can define a 'divergent 
part' of the vacuum energy which is, as known, not uniquely defined. It depends on the kind of 
regularization one has to introduce. For instance, in zeta functional regularization, the regularized 
vacuum energy reads 

^vac(5) = ^^C.^2^ (8) 
n 

where /x is an arbitrary parameter with the dimension of a mass and with a frequency damping 
function it is, 

^^vac((5) = i^u;„e-^-", (9) 

n 

where lu^ are the frequencies of the quantum fluctuations of the electromagnetic field. In our 
problem the spectrum is continuous, but for the moment it is more instructive to keep the notations 
of a discrete spectrum. In zeta functional regularization, the divergent part reads 

i^^c(«) = -^(^+lnM^), (10) 
where we used the notations of [8] in which the heat kernel expansion reads 

^^^^'^ (i^("° + "^^ + ''^* + ---) ■ ^^^^ 
In the scheme with the frequency damping we have 

pdiv/rN 3ao 1 , "1/2 1 ai 1 a2 , . . . 

^vac(<5) - + + + (12) 

The regularizations are removed by s — > resp. 5 — > 0. These formulas follow, for example, from 
section 3.4 in [8] for to = 0. 



4 



The idea of the rcnormahzation is to have in the classical energy -Eciass parameters which can be 
changed in a way to absorb E^^. In the considered model such parameters are the mass m of the 
shell, the frequency cub of the breathing mode, the radius at rest Ra, and the energy £^rest- Now, 
whether this is possible, is a matter of the dependence of the heat kernel coefficients, especially of 
02, on the radius R which is the dynamical variable of the classical system. In the considered, very 
simple model we have only a polynomial dependence on R up to second order in (6). Since we 
assumed adiabaticity for the motion of the shell we do not have a time dependence in a2 so that 
it cannot contain R. Hence, the kinetic energy remains unchanged and, together with it, the mass 
m. Only the remaining parameters, ujb, Rq and -Ercst can be used to accommodate the divergent 
part. In fact, this turns out to be sufficient for the considered model. As it will be seen below, the 
heat kernel coefficients ao, . . . , 02 which enter the divergent part, depend on the radius polynomial 
and at most quadratically. In this way this model is renormalizable. 

It should be mentioned that this scheme is equivalent to the corresponding one in quantum field 
theory with E^!^ in place of the counterterms. Also the interpretation of the renormalization is 
similar. Namely, we argue that the vacuum energy in fact cannot be switched off and what we 
observe are parameters like, for example in QED, electron mass and charge, after renormalization. 

Within this scheme of renormalization, the specific form of the heat kernel coefficients is in- 
significant. The only what one has to bother of is its dependence on R to fit into the freedom 
of redefining the parameters in Sdass- If this is the case, one may define a renormalized vacuum 
energy by means of 

ElZ = lim (Svac(s) - E^Zis)) (13) 
(and the same with 6 in place of s) and one has now to consider 

■E'tot = -£^class + ■E'vac (1'^) 

in place of (7). In this way, the question on how to remove the ultraviolet divergences is answered. 

It remains, however, the question about the uniqueness of his procedure which comes in from 
the parameter fj, in the zeta functional scheme or from the possibility of a redefinition S ^ cS in 
the other scheme. 

In the case of QED at this place one imposes conditions on the analog of -E^™ in a way, the the 
mass and the charge take the values one observes experimentally. 

In our case a similar scheme is conceivable too. A different scheme, suggested in [7], using the 
large mass expansion to fix the ambiguity does not work here since the electromagnetic field is 
masslcss. A way out could be to look for a minimum of the total energy, -Etot, (7), which however 
would imply to take the model (6) seriously. This is not the aim of the present paper. Instead, as 
a normalization condition we demand that in the limit of the plasma frequency f2 — > oo, where the 
matching conditions (2) and (5) turn into that of an ideal conductor, we shall recover the vacuum 
energy of a conducting spherical shell, i.e., just the quantity which was first calculated by Boyer 
in [11]. Indeed, as we will see in the next section, this is possible using the freedom of a finite 
renormalization. 



III. THE JOST FUNCTIONS AND THE HEAT KERNEL COEFFICIENTS FOR THE 

SPHERICAL SHELL 

The electromagnetic field interacting with the plasma shell is defined in the whole space and it has 
a continuous spectrum. In that case the vacuum energy, after the subtraction of the contribution 
of the empty space, can be represented in the form (see Eq.(3.43) in [8]) 

oc ^ 

i?o(.) =-^/.2^^ J dkk'-'^^lnfiiik) (15) 

1 = 1 g 

with v = I + 1/2. The arbitrary parameter /z has the dimension of a mass and fi{k) is the Jost 
function of the corresponding scattering problem. It is determined by the wave equation and the 

matching conditions (2) and (5) respectively for the TE and TM modes. Here we have to consider 
the regular scattering solution which is defined as that solution which for /c ^ turns into the free 
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solution. For r ^ oo, it describes a superposition of incoming and outgoing spherical waves and 

in our model it can be written in the form 

r) = ji{kr)&{R - r) + ^ (/i(A;)/ip^(fcr) + f:{k)h'i'\kr)) e(r - R), (16) 

where j; (a;) = i/7r/2xJ;_|_i/2(a;) and hi^'^\x) = \l tt /2xHj^y^ (x) are the spherical Bessel functions 
and fi{k) and ,fi{k) arc the Jost function and its complex conjugate. For r ^ R these are solutions 
of the radial wave equation. Imposing the matching conditions (2) and (5) on (16), the Jost 
functions can be determined separately for each polarization, 

/i™(fc) ^ l-inkR^ji{kR)hl^\kR), 

/™(fc) = l+z^jl{kR)hl'^'{kR). (17) 
The corresponding formulas for imaginary argument read 

/,™(ifc) = l + jSi{kR)ei{kR), 

f™{zk) = fc2 (^1 - ^ s[{kR)e\{kR)^ , (18) 
where we used the modified Riccati-Bessel functions 

si{x) = \I^Ii+i/2{x), ei{x) = Ki+ii2{x). (19) 

In (17) and (18) we made use of the freedom to multiply the Jost functions by a constant which 
does not influence the vacuum energy (15). 

In zetafunctional regularization, the ultraviolet divergences manifest themselves as poles in the 
of the regularized energy Eq[s), (15). In our case the pole structure reads 



and included the double poles which will appear below in the TM mode. Fortunately, the double 
poles start from A; = 5 and do not influence the renormalization. 

In order to find the poles one has to construct the analytic continuation of i?o (s) into the region 
where the sum and the integral in representation (15) do not converge. For this one may use the 
uniform asymptotic expansion ff^ (ik) of the Jost function for large both, u and k with z = ^ fixed. 
We define 

oo 

E^is) = ^ y dkk'-'^-^lnfriik), (21) 

1 = 1 g 

whose pole contributions coincide with that of Eo{s), (15). In the following subsections we obtain 
the heat kernel coefficients separately for the TE and TM modes. As for the TE modes the 
procedure is well known. One simply inserts the uniform asymptotic expansions of the Bessel 
functions entering (18) and the analytic continuation is an easy task. For the TM mode, however, 
this does not work and a more sophisticated treatment is in order. 



A. The asymptotic expansion for the TE mode 

Directly inserting the known uniform asymptotic expansions of the Bessel modified functions [1] 
into (18) one obtains with k = uz 
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The cj^ are polynomials mt=l /Vl + z"^ (see [3] , Appendix A) and we used 

c™ = ^ (1 _ + 5^4^ ^ ^TE ^ ^'(1- ^') (27 - 553^2 + IGlTi^ - llSSi*^) . (23) 
8 128 

Using this expansion we define the asymptotic part of the logarithm of the Jost function, 

ln/^■"^(^fc)=E§' (24) 



1=1 



with 



Inserting this into (21) we define Eq^'^^(s) 



1=1 Q 

This expression is in a suitable form for the analytic continuation because it can be immediately 
expressed in terms of known functions. In (26), the sum over v results in Hurwitz zeta functions, 

fL'^-' = CH(s;l) (27) 



for the scalar field with the s-wave included and, without the s-wave, 



E^"' = C«(^5|) (28) 



1=1 

for the electromagnetic field. The integration over z can be carried out using 



f 

Jo 



From Eq. (20) we get then the heat kernel coefficients for the TE mode. These are shown in Table I. 
We remind that for the TE polarization there are no double poles which can be confirmed by direct 
inspection in the above formulas. This holds also if higher orders of the asymptotic expansion from 
(22) are included. For the coefficients for the scalar field, i.e., including the s-wave, we confirm 
the results found in [7, 17]. The coefficients for the electromagnetic field become different starting 
from fc = 2. 



k 


/ = 0,1,... 


Z = l,2,... 











1/2 








1 






3/2 






2 







TABLE I: The first few heat kernel coefficients a^. for the TE polarization, including the s-wave in the 
second column and without it in the third column. 
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B. The cisymptotic expansion for the TM mode 



As before in the preceding subsection we insert the asymptotic expansions of the Bessel functions 
into the Jost function (18) and obtain 

The coefficients cj^ are also polynomials in t and we need only 



t 



2 



C 



^^ = -!-(l-et^ + 7t^). (31) 



The direct insertion of this asymptotic expansion into In fi{ik) would produce powers of z in 
the denominator and a term-by-term integration over z would be impossible. Therefore we define 

P = ^^ + ^ (32) 
and perform the asymptotic expansion of the logarithm for large u with fixed p, 

l-/rW-lnp + ln|l + ^|:|-l, (33) 



and define a part of this expansion as 



ln/™'-(^fc)=lnp+^^i^. (34) 

In fact, this is a partial re-summation of the expansion one would obtain acting in the same way 
as in the preceding subsection. Obviously, the pole contributions to E™{s) can be obtained from 
this expansion too. The advantage of (34) is besides allowing for a term-by-term integration that it 
turns into the asymptotic expansion for Neumann boundary conditions in the formal limit i7 oo. 

Inserting now (34) into (21) we obtain for the asymptotic part of the vacuum energy of the TM 
mode the expression 

E^'-is) = M^^EEEEn:n-^-"-%% (^) (35) 

/>! fe=On=l r=0 ^ ' 

with 

oo 

, , COS^S /■ , z2-2-(l + z2)-'=/2 

J (z^ + aVl + z^) 
and the non-zero coefficients Y^^ are given in Eqs. (37): 



(37) 



with Q = QR. As a result of the re-summations these integrals are more complicated than that 
which appeared for the TE modes. 

In order to perform in (35) the analytic continuation in s we would like to expand the integrals 
(36) in a series in powers of a. The direct expansion in the integrand is impossible because of 





= 1, 


= 


Q 

4 ' 








_ Q 

32 ' 


n,! = 


9Q 
16 ' 


Hi 


_ 350 
32 




_ Q 

16' 


y3 _ 

■^3,2 — 


3Q 
8 ' 


n.2 


_ 7Q 
16 


Y4 
^2,2 


_ 

64 ' 


v4 _ 
44 — 


3Q2 
32 ' 


y4 
^6,2 


70^ 

64 
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the behavior for small z. Therefore the idea is to move the integration contour away from passing 

through = 0. 

Splitting the cosine into two exponentials we change the integration in the second part by the 
substitution z — > —z. Both part can be united into one integral over the whole axis and (36) 
becomes 



27r 



v2-2s. 



dz ■ 



1 



(^2 + ayrr^) 



(38) 



Next we move the integration path into the upper half plane. There is a pole of n-th order in 
izq. We move the integration path across this pole which gives an additional contribution which 
results in the second terms in (39), (40) and (41) below. Moving the path further upwards we hit 
the branch cut starting from z = i which originates from vl+^ in (38). To handle the singular 
behavior in z = i we divide the contour into two parts, one is a circle around z = i with a small 
radius e and the other is the path closed to the two branches of the cut with z = ix^ x =\-\-e . . .co. 
Integrating in this integral by parts fc — 2 times, the surface terms just cancel the divergent terms 
coming from the circle around z — i and in the limit e — > we obtain an integral (first terms) and 
an explicit contribution (last terms) in Eqs. (39) and (40), 



and 



An{2k 

•2s- 1 



iZo)"z2-2.i2fc+l 



(n- 1)! 

'{2k-iy. 



(^2 + aV22TT) 



(n-l) 



(i + Z^) 



2^2-2s 



((« + ^2)2 + az^SiTI)" (2i + z2)'=+i/2 



(2fe-l) 



0=0 



'■2k,n 



(a) 



(-2)'= 



47r(2fc - 3) 



/OO 
X) 



+ 



j2«-l 

"^(2fc-2)! 
where we introduced the notations 



i ^ \n_2 — 2sj-2fe 

[Z — IZq ) Z t 



(n-l) 



(^2 + aVi^Tiy 

(i + z2)2-2« 



((i + ^2)2 + azV2iT^)^ {2i + ^;2)fe 



(2fe-2) 



z=0 



(39) 



(40) 



zo = 
Fix) = 

*(a;) = 



ct a j—z 



1 



{—X — ia\/x — ly 
1 



{—X + ia.\/x — 1)^ 



i\/x — 1 



{—X — ia\/x — 1)" {—X + ia\/x — 1)" 



These formulas hold for fc > whereby for = the last term in ^ (a) must be dropped. As 
an example we note for the simplest case with A; = and n = 1, 



2^ 



r 



1 2^ 
+ 



2Q! 



dxx ^~^2,/x'- 

x^+a^{x-l) ' y/^PT4: 



(^-a + Va2 + 4 j 



(41) 



The merit of the representations (39) and (40) is that these can be directly expanded into powers 
of a. For instance, from (41) we get 



-'0,1 
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r(s + i) 



Sa^ r (s + 1) 



(42) 



Now we insert these expansions into '^^(s), (35) with a = VlRl{2v). There the sum over 
the orbital momentum v = I + 1/2 dehvers directly Hurwitz zeta functions, Ci/(a, 3/2), with 
corresponding a. In case the s-wave is included, the result would be expressed in terms C,h{o-, 1/2). 
Keeping the necessary number of contributions we come to 



4(Mi?)2^ r Q{s-\)T{s) 



R 



(43) 



+ 



g(s - i) (302 _ (2s + 3)(2s + 5)((2.s - l)(7s + 'f) + 27)) T (s + 1) 



Ch{s 



i92v^r(.' 



3 3 Qjs-^ 
2'2''^ 4 



Ch(2s + 1, -) + ... 



1 3. 



32 



13, 
Ch{s+-,-) 



+ 



1)2 



l)Q2 + 24(5+f)). 3 3 

-^^i Ch(^+2'2) + --- 



Prom this representation and together with Eq. (20) we obtain the heat kernel coefficients for the 
TM polarization. These are shown in Table II. 

These coefficients can be compared with that obtained for a plane plasma sheet by dividing by 
the area of the sphere, A-kR'^, and taking the limit R^ oo. In fact, these coincide for A; > 1 with 
that obtained in [9]. 



TM 


1 = 0,1,... 


Z = l,2,... 











1/2 




Stt^^^R^ 


1 




-^nR^ 


3/2 


14^3/2 


_10_3/2 


2 




-47rn + 2|n3i?2 



TABLE II: The heat kernel coefficients for TM polarization. We represent the calculation with s-wave 
(second column) and without s-wave (third column). The difference appears starting from k = 3/2. 



From Eqs. (43) together with (20) also the double poles mentioned in the introduction follows. 
It comes at s = — 1/2 from the zeta fimction (s + |: |) in the last line in (35) and the gamma 
function in the left hand side of Eq. (20) and corresponds to heat kernel coefficient 05/2- There is 
no double pole at point s = 1/2 which corresponds to a'^j^ due to factor (s — 1/2) at Ci?(s + 5; |) 
in (43). Therefore the logarithmic contributions for heat kernel expansion starts from flg 



IV. THE RENORMALIZED VACUUM ENERGY 



The renormalization of the vacuum energy is given by Eqs. (13) and (10) with the heat kernel 
coefficients calculated in the preceding section. As discussed in Sec. II, it remains to dispose of the 
freedom of a finite renormalization. It turns out that it is possible to join this with the behavior 
for large f2, i.e. with the limit of the plasma sphere to become an ideal conductor. 

We start from dividing the regularized vacuum energy into two parts, 

£vac(s)=KaT + i^^Ic(«), (44) 

where E^^^ is defined by Eq.(21) and the 'numerical' parts of the energy are defined by 

^vaT = -^E^/ dkk^{\nm)-\nfr{ik)). (45) 
;=i Q 



10 



In this way the asymptotic part of the logarithm of the Jost function, \nfl^^{ik), is subtracted in 
i?"ac" ^-i^d added back in iJ^ac- general, the definition of ln/f"(ifc) is not unique. The only one 
has to ensure in any case is that the sum and the integral in (45) do converge if one puts s = 
there. We did this in Eq.(45). This part is called 'numerical' since it allows for a direct numerical 
evaluation. 

In the preceding section we used the expressions (24) and (34) for In f^{ik) for the calculation 
of the heat kernel coefficients. In this section we take the same expression for the TM mode (34) 
and different expression for TE mode, Eq.(46) below. In In/, '"""(ifc) we make a re-sumniation 
like that in In f™''^{ik), (34), in order to archive also for \nf^^''^{ik) the property in the formal 
limit ^ oo to turn into the corresponding expression for an ideal conductor. 

With this motivation we define 

ln/™-(^fc)=ln«,+ ^^, (46) 

where 

«, = ! + —. (47) 

We would like to stress again that a redefinition of In /,^^ {ik) is merely a redistribution of contribu- 
tions between and E^l.{s). With the definitions (46) and (34) we archived that in the 
limit O — > c» must turn into the corresponding ideal conductor expressions. This follows because 
in (45) both, the Jost function and the part (46) of its asymptotic expansion, do that and because 
in addition the integral and the sum are convergent. Indeed, from the numerical evaluation we got 

lim E--'-- = 

— s-CSO It 

which is the same as if one takes conductor boimdary conditions from the beginning. With the 
other part, E^^{s), this is not such simple. The corresponding calculations are carried out in the 
Appendix and the result is 



E^TE,as/„^ _ "2 
^vac («) - 

„T 

i;™'-(s) = 



TE 



327r2 



1 n 

- - 2 In — 
s 2/i 

1 ^, ^ 

- - 2 In — 
s 2ii 



+ 1800^^36^+^-^ 



These expressions are sums of a divergent part (it is proportional to the heat kernel coefficient 02 
as expected), two terms growing with fl and an 'analytical' part, 



4 



1=1 ;=i 

=2 ;=i 



All quantities entering -B™' ^" ^^id defined in the Appendix. The analytical parts have 

the 'necessary' limit for 00, 

lim £™,an ^ J]_ 

Q^OO ™^ 128i?' 

Jim i^™-- = (51) 



With Eqs. (49) and the property (51) we have all information we need to complete the rcnor- 
malization. We remind the discussion in Sec. II that all terms which are proportional to R, B? 
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or which do not depend on R can be removed by a redefinition of the parameters in the classical 
part. In (49) this concerns all except the last ones, i?™'™ and E™''^'^. That means, that we not 
only can remove the contribution proportional to 02, but also the terms growing with fl. For this 
reasons we define the renormalized vacuum energies by 



jpTE, ron 

pTM, ren 
vac 



pTE, num , pTE, an 
-^vac ~i~ -'^vac ' 

pTM,num , pTM, an 
vac ~r vac 



(52) 



With these formulas we completed the model consisting of the classical energy and the vacuum 
energy which is the sum of the two contributions in (52). The main merit of this vacuum energy is 
that it turns for 00 into the ideal conductor limit. Using the formulas for the Jost functions 
and their asymptotic parts and also the formulas in the Appendix, it is possible to evaluate this 
vacuum energy numerically. The results are shown in the figures IV and IV for the dimensionless 
function £ defined by 



-'^vac 



as functions of their arguments x = flR. 



£{nR) 

R 



(53) 



£(n«) 



£(nff) 



1" -0.005 



0.2/ 0.4 0.6 0. 



OR 



FIG. 1: The function £{Q,R) = RE^^ plotted as function of Q,. For large O it tends to the ideal conductor 
limit, limfiH^oo — 0.0046 (left panel). For small Q. (right panel) it takes negative values and decreases as 
£{ilR) 0.0589VI1R. 




FIG. 2: The renormalized vacuum energy divided by O, iSvaS/f^, as a function of R. For large R it 
approaches the ideal conductor limit and for small radii, R < it becomes attractive. 

It is also interesting to consider the limit of small argument of these functions which is equivalent 
to a small plasma frequency O. For the TE mode the main contributions come from the integrals 
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An 



For the TM mode we get accordingly 

j2(nR) = ^^^^^^^}^^^VnR+ OiQRlnflR), 
647r2v2 

JsiQR) = ^fiR ^^" ^^5^(2) +o{QRlnnR), 

8v 2 

Ji{D.R) = --^^ +0{nRlnnR), 

~ \/nR 

JaiflR) = +0(0i? In f2i?), 

8 

where Cni^') is the Ricmann zeta function. The numerical parts are ~ 0(f2) which is easy to show 
using the next term of asymptotic expansion of the Jost functions. In this way we obtain the 
following behavior of the vacuum energy for small O, 

pTM.ren, ^ / 3(4-V2)C(|) (l-V2)C(i) 1^ 

= -0.0598W^. (54) 
V R 

Therefore the main contribution for the energy comes from the TM polarization. The energy tends 
to zero proportional to a/TI. The same behavior was observed in Refs. [9] and [5]. 



V. CONCLUSION 



In the foregoing sections we considered the vacuum energy of the electromagnetic field interacting 
with a spherical plasma shell. We calculated the heat kernel coefficients for both polarizations. 
For the TE case the standard methods apply, for the TM case a re-expansion of the asymptotic 
expansion of the logarithm of the Jost function was helpful. It turned out that the vacuum energy 
in zeta functional regularization and, with it the corresponding zeta function, have double poles. 
This implies that the corresponding spectral problem is not elliptic. On the other hand, at least 
based on the calculations carried out in this paper, there is nothing which would diminish the 
reasonability of this model. 

A basic concern of this paper is to construct a model allowing for a physically meaningful 
interpretation of the renormalization. We considered with the breathing mode of the shell the 
simplest model for the classical motion of the shell. It turned out that this model is able to 
accommodate all renormalizations which we were like to carry out. These are the removal of the 
pole in s, i.e., of the ultraviolet divergence, and the removal of all contributions growing together 
with the plasma frequency f2. It should be mentioned that this includes also the removal of the 
arbitrary constant which came in with the regularization. The nontrivial statement which allowed 
for doing so is that the dependence on the radius R of all these contributions is polynomial not 
exceeding R^. 

It would be interesting to investigate the question whether this procedure can be carried out 
also for more general deformations of the shell. In principle, most ingredients for such a calculation 
are available. Especially, the heat kernel coefficients for the TE modes can be taken from [17]. It 
would remain to calculate the coefficients for the TM modes. 

Concerning the arbitrariness of the normalization procedure we would like to mention that 
the removal of the contributions growing together with Q can be considered as a normalization 
condition. It ensures the uniqueness of the renormalized vacuum energy and makes this model 
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physically meaningful. In this way. the gap between the renormalization procedure in quantum 
field theory in smooth background fields and the removal of divergences of the Casimir energy in 
the background of boundaries, as suggested for example in [4] (section 6.5), is narrowed. At once 
in this way the much discussed vacuum energy of a conducting spherical shell now appears as a 
limiting case of a slightly more physical model. 
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VI. APPENDIX 

In this appendix we perform the analytic continuation in the asymptotic parts of the vacuum 
energy defined in section IV. For that, we use the following integral representations of some sums, 

_ 7ra3-2^ dyy^-^' -acosTTs + ysinns 3 



1=0 

,2-2s ^„2-2s foo , ,„,2-2s 



1=0 

Eu'^~^^ Tra^"^** „ /■°° -acosTTS -I- wsinTTS 1 ^ 
= r^; 2 / ; , - <Us < 1, 
i_^t^ + a sin27rs Jo + a? l + e^'^v 2 

^ 27r(s-i)a-2^ r"^ dyy-^' 2ay sm.-K s - {a? - y'^) cos -ks ^ 



y 

dyy^~'^'^ a simrs + y COS TTS c0 3 



fr^i^ + a)^ sin 27rs Jo iv^ + a'^Y 1 + e^''^ 



1 r dy I 

wa 



1=0 + "^^^ -^0 + "'^ ^^^^^ ' 

which were obtained using the Abel-Plana formula in the form 

f{-iy) - fi-iy) 



1 POO pc 

1=0 ^ -^0 



dx 



1 + e^'^y 



A. TE polarization, h = Qt/2 

The expression for asymptotic part ii^™' appearing from inserting (46) into (26), after carrying 
out the differentiation, has the following form 



OO POO 



f Qi3 (3 - 30^2 + 35*4) g2^6 ^ _ 6^2 ^ 5^4^ I 

with w = 1 + Qtjlv and Q = QR. We perform the calculations separately for the contribution 
from each power of u using the formulas (55) given above. 

r -It . '^OSTTS p^2s ,,2-2s f°° j 2-2s Q^^ 



nOO 



2wi' 



14 



,2-2s 



1=1 



1 



11 



487r 

cos TTS 
TT 



1 ^, 

- - 2 In — 
s 2/x 



0^ 



+ ^ + Vi(g) + V2(g) + Vi(Q), 



IGtt 



/•OO 



-2s 



+ 6 



967r 



1 ^, 

- - 2 In — 

s 2fj, 



Q 



50407r 



+ V3(Q)+V2(Q), 



M : -'^{,,Rf^Y,y^-^^ dzz 



32«;2 



3Q 
5607r 



-V4(Q)+V3(Q), 



where the following integrals were introduced, 



Vi 

V2 
V3 
V4 

Vi(Q) 



271- io 1 + e^'^?' 



In 



1 + 



4y_ 



.21 



1 + e^^y Jo 



1 + e'^^'V y2 _^ Q2i2/4 l+xt' 

ydy 



Q_ 

St: Jo 1 + e^^v Jo 

Q_ y^dy I 

16 Jo cosh^ Try Jo 



y2 + Q2i2/4 ' 
~ ^^2^5(1 _ 6t2 + 5f4)(ix 



y2 + Q2i2/4 



47r 



i2 /-oo ^2+4 



l + Qt 477 



(l-ln2Q) 



( 7r - 4Q + i7rV-l + 4Q2 + 2^-1 + 4Q2arctanh 



2g 



V-1 + 4Q2_ 



V2(Q) = / da;a;2t 



3 - 30i2 + 2,bt^ 
l + Qt ' 



^^^^^ = s^rX '""^ (1 + Q^)2 • 



We note the following expressions which are necessary to consider the ideal conductor limit, 



jToo^^(^) = 8' 



All other integrals vanish in this limit. 
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B. TM polarization, a = Q/2tz^ 

The expression for asymptotic part E™' appearing from inserting (34) into (26), after carrying 
out the differentiation, has the following form 



iTM, as 



2 cos ITS 

'l Qt_ Qt ( (35t^ - 18^^ + 1) e lAt^ - 12e + 2 \ Q'^ {7t^ - 6t^ + l) 
p 4pi/ 32u^ \ P I 



p = z"^ + Q/2yt. As before, we perform the calculations separately for the contribution from each 
power of V using the formulas (55) given above: 



^ h Jo P 
-^(,i?)- rdzz-^^y^'^ 



HQ 



2 cos TTS 
TT 

, COS TTS 



(/Zi?) 



2s 



1 r., ^ 

- - 2 In — 
s 2/i 



dzz 



IOOtt 

2-2s 



+ MQ) + J2{Q) + MQ), 



Avp 



= _Q^(^ij)2. rdzz--Hy^^^ 

27r Wo tit ^ + " 



iig 

487r 967r 

2 cos TTS 



1 

- - 2 In — 

s 2/x 



+ ^-\mQ) + Js{Q) + J2{Q), 



-Q 



Jo 



, Qt 

32i/3 



^^(35^^ - ISt^ + 1) 2(7t^-6t^ + l) 



COS TTS 



23Q 
96^ 

2 COS TTS 



1 , fl 

In — 

s 2/i 



dzz^-'H 
547(3 



t^(35t^- 18^^ + 1) ^ u-"^' 2(7t^ - 6^^ + 1) 2^- 



l-2s 



7207r 



+ J4((3) + J5(Q) + J3(<3), 



[W : -i^^(^i?)2«^^2-2.^ ^^^2 



^Q2^2(7^4_g^2^1^ 



64z/4p2 



r dzz-H\n^ -&e + i)y 



HQ 

'2407r 



+ J6(Q) + J4(Q), 



where the following integrals were introduced. 



:^i(0) 

J2{Q) 

MQ) 
MQ) 



Q 



3 (-oo 



zdz 



arctanhVl — z 

, h In — 

2 

x^tdx 1 



16n Jq 1 + e'^Q^/^ 
2Q /"^ 



/■°° y^dy 1-°° _ 



x^t'^dx 



= ^ 

t: Jo 1 + e^^y Jo a;4y2 + ^2/4^2 i + xt' 

= ^ r dxt^x\35t'' - ISt^ + 1) /' 
Stt Jo Jo 



ydy 
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MQ) = f / - 6<V 1) / - 

JeiQ) = ^ r dxt'x^n^-Qe + l) I 





MQ) = 

MQ) = 

MQ) = 

MQ) = 



xtdx 



2t: Jo {x"^ +Q/t){x + Q)l + xt' 
^2 1-°^ dx Q , „^ 



Q2 /-o 
47rio 



Stt 



dxx 



x'^+Q/t 4n 

ri3(35t4 _ 18^2 _^ 2t(7t'* - Gt^ + 1) 



Stt 



L 



x2 + g/t 

2 '■~d:ct2x2(7f4- 6^2 + 1) 



(x2 + g/t)2 



(x2 + g/t)2 

We note the following expressfons which arc necessary to consider the ideal conductor limit, 

lim MQ) = 



lim MQ) 

Q— >oo 



7 

256' 



All other integrals vanish in this limit. 
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